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We show that certain little Higgs models with symmetry breaking SU(N) — » SO(N) for N > 4 
admit topologically stable solitons that may contribute to cosmological dark matter. We have 
constructed a spherically symmetric soliton and estimated its mass in the case of 5(7(5) — > SO(5). 
Its lower bound is found to be around 10.3 TeV. Whether this particle is a fermion or a boson 
depends on the value of an integer-valued parameter of the underlying theory, analogous to the 
number of colors of QCD. In either case, the particle is neutral. If it is a fermion, it is a Majorana 
particle, which could take part in a seesaw mechanism for neutrino masses. 



The hot Big Bang model of cosmology tells us that the 
Universe cooled down from a primordial hot and dense 
state to the present state of galaxies and other large- 
scale structures with a mean temperature of about 2.73 
K [T]. Certain theories of grand unification predict that 
the Universe underwent through a series of phase tran- 
sitions as it cooled down, similar to what we observe in 
condensed matter systems. Phase transitions in the early 
Universe can give rise to certain stable configurations of 
matter known as topological defects. Different types of 
topological defects can arise depending on the symmetry 
breaking mechanism of the underlying field theory [2\. 
They can appear in the forms of magnetic monopoles, 
cosmic strings, domain walls, textures and skyrmions [S]. 

A field theory described by a continuous symmetry 
group G, when spontaneously breaks down to a subgroup 
H C G, the space of all accessible vacua of the theory 
called the vacuum manifold is defined to be the space 
of cosets of H in G. The theory possesses a topolog- 
ical defect if some homotopy group of the coset space 
ir d (M = G/H) is nontrivial. When d = 0,1,2 the de- 
fects respectively are domain walls, strings (vortices) and 
magnetic monopoles or textures. The case d = 3, which 
plays a major role in this paper, gives rise to point-like 
topological defects called skyrmions. 

Recently there has been much interst in a class of field 
theoretic models called the Little Higgs models [H [5j |6] 
in the context of weak scale symmetry breaking. These 
models provide a new logical possibility for natural elec- 
troweak symmetry breaking and a new partial resolution 
of the hierarchy problem in elementary particle physics. 
Introduction of new symmetries at the TeV scale by these 
models provides the cancellation of all quadratically di- 
vergent contributions to the Higgs mass at the one-loop 
level and pushes up the hierarchy problem to an energy 
scale of around 10 TeV. Little Higgs models have gen- 
erated a lot of interest since any potential candidate to 
solve the hierarchy problem deserves serious attention. 

Among the many possible ways of implementing the 
little Higgs paradigm, the littlest Higgs model [5] is the 
simplest and most economical. This theory introduces a 
weakly coupled new physics at TeV energies, stabilizes 



the electroweak scale with a naturally light Higgs and is 
the smallest extension of standard model to date. 

In this paper we address the interesting new possibil- 
ity of bridging the natural electroweak symmetry break- 
ing and cosmological dark matter - the non-baryonic, 
non-relativistic and weakly interacting matter that con- 
stitutes about 22% of matter in the Universe. Since 
the Higgs particles appear as pseudo Nambu-Goldstone 
bosons in little Higgs models, skyrmion solutions that 
are stable and electrically neutral can also come out quite 
generically. We demonstrate the existance of a point-like, 
electrically neutral and topologically stable structure; a 
particle with a Zi charge; which could be a viable dark 
matter candidate. Its mass, estimated in the context of 
littlest Higgs model with T parity 7 , is found to have a 
lower bound of around 10.3 TeV which is well below the 
unitarity limit [5] of viable dark matter particles. (Exis- 
tence of other topological defects in the little Higgs model 
was investigated in [9].) 

We start with a class of non-renormalizable effective 
field theories for pseudo Nambu-Goldstone bosons, in 
which a symmetry group SU(N) is broken down to its 
real subgroup SO(N) for N > 4. The case N = 5 is of 
most interest, as it appears in the little Higgs models. 

At small energy scale compared to the symmetry 
breaking scale, the effective action has the form 

Si = y J d 4 xtr d^d^ + ■■■ , (1) 

where / is a parameter with dimension of energy and 
$ is a scalar field given by a differentiable map from 
the Minkowsky space R 1,3 to a nonlinear target manifold 
M N , 

$ : R 1 * 3 -> M N . (2) 

The target manifold A4n is the subset of symmetric ma- 
trices in SU(N) 

M N = {$|$ = $ T , = 1, det$ = 1}. (3) 

It has a global symmetry $ — * g$g T , g £ SU(N). Any 
$ £ M. n can be reduced to the identity by this transfor- 
mation [10]. That is, there is a g £ SU(N) such that 
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<f> = gg T . If we change g >— > <?/i, with /i S SO(N), 
the product g<7 T is unchanged. Thus we can identify 
M N = SU(N)/SO(N). The canonical projection to the 
cosets is p : SU(N) — * Mjv> = ff5 T - 

At spatial infinity, the field <I>(:e) must approach a con- 
stant; the choice of this constant among all matrices sat- 
isfying Eq.(|3| will break the symmetry SU(N) down to 
its real subgroup SO(N). The parameter / sets the scale 
of the symmetry breaking; in the little Higgs models it 
is expected to be a of the order of a TeV. The dots in 
Eq.([T|) indicate that we are ignoring higher derivative 
terms, which are expected to be unimportant in the limit 
of 'low' energies; that is, energies of the order of /. 

Among the higher derivative terms we can add a new 
term that does not change the hyperbolic nature of the 
field equations and is still second order in time. This is 
the "Skyrme term" [3] given by 



32e z J 
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The value of the dimensionless constant e depends on 
the details of the renormalizable theory of which Eq. 
is the effective action. We will see that in the presence 
of this term, the effective action supports a topological 
soliton, whose mass is proportional to M = fl/e, with I 



given in Eq. ( 22 ) 



The more familiar Skyrme model 112] is for the 
spontaneous breakdown of the symmetry SU(2) x SU (2) 
to SU(2). The Nambu-Goldstone bosons are then the 
7r-mesons. The action of the Skyrme model is then 



f 2 

■I 7T 



d 4 xtr d^gd^g^ 



+ i / d4xtT i dfl 9,d v g^Wg,d v g^ + ■ ■ ■ (5) 

Closer in spirit to these papers are the Hopf soliton [T3] 
(the case of A4 2 ) and even more so, the model studied in 
[T4"] . which is the case of Ai 3 . 



I. TOPOLOGICAL CONSERVED CHARGE 

A continuous function $ : i? 3 — > A4n that approaches 
a constant at infinity can also be thought of as a map 
$ : S" 3 — > A4n by identifying the points at infinity. The 
homotopy group tt 3 (M.n) has as elements the equivalence 
classes of such maps that can be deformed continuously 
into each other. It is well-known that [15] 

t: 3 (M 2 ) = Z,tt 3 (M 3 ) = Z 4 ,tt 3 (M n ) = Z 2 ,N> 4. (6) 

The case N = 3 was studied in a different context some 
years ago [14]. We will focus here on the case N > 4, 
which includes the little Higgs models. There is just one 
non-trivial equivalence class of maps $ : R 3 — ► Mn', 
we will need to determine which representative of this 
class has the least energy. The non-trivial element of 
7r 3 (SO(N) I SO(N)) is just the projection of the genera- 
tor Of 7T 3 (SU(N)). 



II. A SPHERICALLY SYMMETRIC ANSATZ 

Recall Skyrme 's spherically symmetric ("hedgehog") 
ansatz for a soliton of winding number one: 

92(x) = e "-M0 w ( ) = w (oo) = Q. (7) 

The boundary conditions on uj ensure that the limits 
at r = 0, 00 are direction independent: 



32(00) = 1 2 , 32(0) 



-Is 



(8) 



This ansatz is spherically symmetric in the sense that 
a rotation in space can be compensated by the adjoint 
action of SU{2): 



fj2 



(R(A)x) = Ag 2 {x)A\ 



(9) 



where R : SU(2) — > SO(3) is the usual homomorphism. 
The obvious topologically non-trivial map g 2 g 2 is n °t 
spherically symmetric: it is just cylindrically symmetric 
around the a; 2 -axis. This is because the representative 
A of the rotation matrix does not cancel (is not orthogo- 
nal) unless the rotation is around the x 2 axis. In fact, the 
energy minimizing configuration in M. 2 is only cylindri- 
cally symmetric |13j . If there is a spherically symmetric 
configuration, it is likely to have less energy. 

There |14j is another spherically symmetric map g 3 : 
R? — > SU(3), which interpolates between the identity at 
infinity and a cube root of unity at the origin: 



93(00) = 1 3 , 53(0) = e 3 1 3 . 



(10) 



It is a generator or tv 3 (SU(3)). To construct it we start 
with the spherically symmetric ansatz 



[93(x)]kl = A(r)[S k i - x k xi}+ B(r)e k i n x r , 
+C(r)x k xi, £i = Xi/\x\, 



(11) 



with the constraints |C| = 1,A*B = B*A, \A\ 2 + \B\ 2 = 1 
to be unitary and C(A 2 + B 2 ) = 1 to have determinant 
one. Under the action g 3 (x) —> Rg 3 (Rx)R T this is spher- 
ically symmetric. 
So we get 

A(r) = e-J^cosaW, C(r) = e lx{r \ 
B(r) = e~5 x(r) sina(r). (12) 
The boundary conditions 

X(oo) = 0, x(0) = 2tt/3, a(oo) = 0, a(0) = tt (13) 

ensure that the winding number is one. Computing 
*a(«) =g 3 {x)gl{x), 

[*s(x)]iu - e-^Skt + [e 2i ^ - e-^xkxt. (14) 

Finally, we can embed in M n to get a spherically sym- 
metric representative for the generator of it 3 (M.n) for 
N > 4: 



* 3 (a:) 
ljv_a 



(15) 
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For $ N (x) = g N (x)gJj(x) with 

9n{x) = I o 1 



N-3 



(16) 



and gw ■ R 3 — > SU(N) has winding number one by the 
above construction. The configuration is spherically sym- 
metric under the action $(2;) i— > RQ(Rx)R T . 



III. MINIMUM ENERGY SOLITON 

The mass of the solition in the theory with action Si - 
S2 will be the minimum of the energy 
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where 



/ 2 (*) 



d 3 xtr di$di&, 



1 

32 



d 3 xtr [a i $,5 j $ t ][5 i $,a 7 -$ 
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(17) 

(18) 
(19) 



among all functions $ : R 3 — > .Mat equivalent to the 
non-trivial element of ^(.Mat). Since this topological 
charge is valued in Z2, the topological soliton and its anti- 
particle are identical. In the absence of other interactions 
a single such soliton will be stable. Their number is not 
conserved - a pair of them can annihilate when they come 
in contact with each other. 

As with skyrmions in QCD, it is clear that under a 
scaling Q\(x) — Q(\x), the two terms in the energy scale 
opposite to each other: 



/i(* A ) = A" Ji($), / 2 ($ A ) = Ala(S). 



(20) 



Minimizing in the scale parameter, we see that the min- 
imum energy will be proportional to f/e: 



H min = (f/e)y/l 1 (Q)I 2 m. (21) 



We can make a variational estimate for the constant 

./ ■ : mil) v /,!<!>) A! <!>>. (22) 



min v / /i($)/ 2 (^)- 



The minimizing configuration should be invariant under 
the simultaneous rotation of the coordinate x and a ro- 
tation by some SU(2) subgroup (analogous to isospin in 
the Skyrme model) of SU(N). 

It is not difficult to make an estimate for the soliton 
mass M. Substituting the spherically symmetric ansatz 
Eq.(15l in Eq.(17) and after some calculation we find 



E(X) 



7T / / 



3r 2 x' 2 +4(l-cos3x) dr 
n I ^(1 -cos3x) 2 dr 
f 2x' 2 

7r / — 2— (3 — cos3x — 2cos6x)dr (23) 




FIG. 1: The solution for x( r ) i n units where e = / = 1. 



We can find the minimum of energy in two ways, (i) 
by taking a variational ansatz for x(r) or (ii) by solving 
-E'(x) numerically. The variational ansatz gives an answer 
almost as good as the numerical solution. 

We tried the following "stereographic" ansatz for x( r ) 



v(r) = — arctan — 
AW 3 V r n J 



(24) 



for n = 1,2,3,4. They satisfy the boundary conditions 
given in Eq.(13l. The lowest value for energy was ob- 

H?. The value of the 

ef 

105 -. The numerical 



taincd for n = 2 and R = R Q ? 

minimum energy is E(Rq) = M 
solution of the differential equation for \ gives a slightly 
lower value of energy close to the variational ansatz: 



E = 102. 



(25) 



We plot the solution in Fig. [T]in units where e = / = 1; 
the dashed curve is the variational ansatz with n — 2 
and the solid curve is the numerical solution. We need 
an estimate for the value of the dimensionless constant e 
as well as / to get a number for the mass of the soliton 
M. Precision electroweak constraints put a lower bound 
on the symmetry breaking scale / typically of about 500 
GcV [16 for little Higgs models with T parity 0. The 
Skyrme constant is in principle determined by the un- 
derlying renormalizable theory of which the little Higgs 
model is the effective theory. At the moment we do not 
know what this effective theory is; even if we knew it, we 
do not yet know how to compute such constants. But it 
is reasonable to expect that e will have the same order of 
magnitude as for QCD; this is the best we can do with 
our current knowledge. In QCD e ~ 5, as we can deduce 
from the value of the nucleon mass. With these values we 
get an estimate (a lower bound for e = 5) for the mass 
of the soliton 



M > 10.3 TeV. 



(26) 



Since the mass of this particle is below the unitarity 
bound (< 340 TeV) 0, it cannot be excluded from the 
list of viable dark matter candidates. Possible cosmo- 
logical implications such as relic abundance, decay |22j 
and annihilation cross-sections of these particles should 
be explored. 
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The coefficient N c of the Wess-Zumino-Witten term 
(which is equal to the number of colors of QCD) deter- 
mines whether the baryon is a boson or a fermion: for 
odd N c it is a fermion and for even N c it is a boson. We 
do not yet know if the analogous parameter in the little 
Higgs models is even or odd: both possibilities would give 
the same effective theory at the electroweak scale. When 
our skyrmion is a boson, it can be represented by a real 
scalar field S whose couplings have the discrete symmetry 
S — > —S. The conserved quantity associated to this sym- 
metry is the topological charge of the little Higgs model. 
The phcnemenological consequences of such scalars have 
been investigated in |17j . When our particle is a fermion, 
it is a Majorana particle. In this case it could be the 
fermion responsible for the neutrino masses in a seesaw 
mechanism |18j . To flesh out this idea, we need to un- 
derstand the mixing matrix of the topolgical soliton with 
neutrinos, induced by the anomalous coupling of neutri- 
nos of the bosons of the little Higgs models [151120] . 

Note added 



While this manuscript was in preparation a preprint with 
a significant overlap with the work presented here ap- 
peared in the arXiv [21] . 
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